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ABSTRACT 

We present an analysis of a pair of Suzaku spectra of the soft X-ray background (SXRB), obtained 
from pointings on and off a nearby shadowing filament in the southern Galactic hemisphere. Because 
of the different Galactic column densities in the two pointing directions, the observed emission from the 
Galactic halo has a different shape in the two spectra. We make use of this difference when modeling 
the spectra to separate the absorbed halo emission from the unabsorbed foreground emission from the 
Local Bubble (LB). The temperatures and emission measures we obtain are significantly different from 
those determined from an earlier analysis of XMM-Newton spectra from the same pointing directions. 
We attribute this difference to the presence of previously unrecognized solar wind charge exchange 
(SWCX) contamination in the XMM-Newton spectra, possibly due to a localized enhancement in the 
solar wind moving across the line of sight. Contemporaneous solar wind data from ACE show nothing 
unusual during the course of the XMM-Newton observations. Our results therefore suggest that simply 
examining contemporaneous solar wind data might be inadequate for determining if a spectrum of 
the SXRB is contaminated by SWCX emission. If our Suzaku spectra are not badly contaminated 
by SWCX emission, our best-fitting LB model gives a temperature of log(TLB/K) = 5.98tg;o4 and a 
pressure of pLu/k = 13,100-16,100 cm^^ K. These values are lower than those obtained from other 
recent observations of the LB, suggesting the LB may not be isothermal and may not be in pressure 
equilibrium. Our halo modeling, meanwhile, suggests that neon may be enhanced relative to oxygen 
and iron, possibly because oxygen and iron are partly in dust. 

Subject headings: Galaxy:halo — Sun: solar wind — X-rays: diffuse background — X-rays: ISM 



1. INTRODUCTION 

The diffuse soft X-ray background (SXRB), which is 
observed in all directions in the ^0.1-2 keV band, is 
composed of emission from several different components. 
For many years, the observed 1/4-keV emission was 
believed to originate from the Local Bubble (LB), 
a cavity in the local interstellar medi um (ISM) of 
~100 pc radius filled with ~10^ K gas (Sanders et aLl 
1977t ICox fc R evnolds [l987l : iMcC ammon & Sander! 
19901 ISnowden et al.l [1990). However, the discovery 



of shadows in the 1/4-keV background with ROSAT 
showed that ^^50% of the 1/4 keV emission originates 
from beyond the LB (|Burrows fc Mendenhalll Il991l : 
ISnowden et al][l99l . This more distant emission origi- 
nates from the Galactic halo, which conta ins hot gas with 
temp eratures log(Tiiaio/K) ~ 6.0 -6.5 (ISnowden et al.l 
1998t iKuntz fc Snow deJ [2OOOI: ISmith et aL 
Galeazzi et 311120071 : iHenlev. Shelton. fc Kuntd [200 



As the halo gas is hotter than the LB gas, it also emits at 
higher energies, up to ~1 keV. Above ^1 keV the X-ray 
background is extragalactic in origin, and is due to un- 
resolv ed active galactic nuclei (AGN; iMushotzkv et al.l 
[2001 . 

X-ray spectroscopy of the SXRB can, in principle, de- 
termine the thermal properties, ionization state, and 
chemical abundances of the hot gas in the Galaxy. 
These properties give clues to the origin of the hot 
gas, which is currently uncertain. However, to de- 
termine the physical properties of the hot Galactic 
gas, one must first decompose the SXRB into its con- 
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stituents. This is achieved using a technique called 
"shadowing" , which makes use of the above-mentioned 
shadows cast in the SXRB by cool clouds of gas be- 
tween the Earth and the halo. Low-spectral-resolution 
ROSAT observations of the SXRB were decomposed 
into their foreground and background components by 
modeling the intensity variation due to the varying 
absorption column density on a nd aroun d shadow- 
ing c lo uds (iBurrows fc Mendenhalll[l9 91: Sno wden et al.l 
flggil, '200a! 'Snowden. McCa mmon. fc Verteil [l99l 
.Kuntz, Snowden, & Verier 199"$! The same tech- 
nique was used to decompose R OSAT All-Sky Surve; 
data over large areas o f the sky (jSnowden et al.l fl99 
IKuntz fc Snowdenll200l . 

With higher resolution spectra, such as those from 
the CCD cameras onboard XMM-Newton or Suzaku, it 
is possible to decompose the SXRB into its foreground 
and background components spectroscopically. This is 
achieved using one spectrum toward a shadowing cloud, 
and one toward a pointing to the side of the cloud. The 
spectral shape of the absorbed background component 
(and hence of the overall spectrum) will differ between 
the two directions, because of the different absorbing col- 
umn densities. Therefore, by fitting a suitable multicom- 
ponent model simultaneously to the two spectra, one 
can separate out the foreground and background emis- 
sion components. Such a model will typically consist of 
an unabsorbed single-temperature (IT) thermal plasma 
model for the LB, an absorbed thermal plasma model for 
the Galactic halo, and an absorbed power-law for the ex- 
tragalactic background. The Galactic halo model could 
be a IT model, a two-temperature {2T) model, or a dif- 
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ferential emission measu re (DEM) model (jGaleazzi et alj 
l2007t iHenlev et al.l[2007l : S. J. Lei et al., in preparation). 

Recent work has shown that there is an additional 
complication, as X-ray emission can originate within the 
solar syst e m, via solar w ind charge exchange (SWCX; 
IC^ [19981 : [C raven iHooO). In this process, highly ion- 
ized species in the solar wind interact with neutral atoms 
within the solar system. An electron transfers from 
a neutral atom into an excited energy level of a solar 
wind ion, which then decays radiatively, emitting an X- 
ray photon. The neutral atoms may be in the outer 
reaches of the Earth's atmosphere (giving rise to geo- 
coronal emission), or they may be in interstellar ma- 
terial flowing through the solar system (giving rise to 
heliospheric emission). It has been estimated that the 
heliospheric emission may contribute up to ^^50% of the 
observed soft X-ray flux (| Craven i [200l . The geocoro- 
nal emission is typically an order of magnitude fainter, 
but during solar wind enhancements it c an be of similar 
brig htness to the heliospheric emission (IWargelin et al.l 
l2004i) . SWCX line emission has been observ ed with 
Chan dra, XMM-Newton, and Su zaku (fWargcli n et al 
2004t ISnowden. Collier, fc Kulrti l2004l : iFujimoto et al 



20071) . As the SWCX emission is time varying, it cannot 



easily be modeled out of a spectrum of the SXRB. 

If SWCX contamination is not taken into account, 
analyses of SXRB spectra will yield incorrect results for 
the LB and halo g as. This pape r conta ins a demonstra- 
tion of this fact. iHenlev et al.l ()2007[ ) analyzed a pair 
of XMM-Newton spectra of the SXRB using the pre- 
viously described shadowing technique. One spectrum 
was from a direction toward a nearby shadowing filament 
in the southern Gala ctic hemisphere [d — 230 ± 30 pc; 
iPenprase et al.l[l998[ ). while the other was from a direc- 
tion ~2° away. The filament and the shadow it casts in 
the 1/4-keV background are shown in Figure [T] Contem- 
poraneous solar wind data from the Advanced Composi- 
tion Explorer [ACE) showed that the solar wind was 
steady during the XMM-Newton observations, without 
any fiares or spikes. The proton flux was slightly lower 
than average, and the oxyg en ion ratios were fa irly typ- 
ical. These observations led IHenlev et al.l (|2007f ) to con- 
clude that their spectra were unlikely to be severely con- 
taminated by SWCX emission. Using a 2T halo model, 
they obtained a LB temperature of log(TLB/K) = 6.06 
and halo temperatures of log(rhaio/K) = 5.93 and 6.43. 
The LB temperature and the hotter halo temperature are 
in good agreement with other recent meas urements of the 
SXRB using XMM-Ne wton and Suzaku ([Galeazzi et al.l 
l2007HSmith et al.ll2007D. and with analysis of the ROSAT 
All-Sky Survey (|Kuntz fc Snowde n 2000). 

We have obtain ed spectra o f the S XRB from the 
same directions as IHenlev et "all 's (|2007f ) XMM-Newton 
spectra with the X -ray Imaging Spectrometer (XIS; 
iKovama et al.l I2007D onboar d the Suzaku X-ray obser- 
vatorv (|Mitsuda et aIll2007D . The XIS is an excellent 
tool for studying the SXRB. due to its low non- X-ray 
background and good spectral resolution. Our Suzaku 
pointing directions are shown in Figure [TJ We analyze 
our Suz aku spectra using th e same shadowing technique 
used bv lHenley et al.l (l2007| ). We find that there is poor 
agreement between the results of our Suzaku analy sis and 
the re sults of the XMM-Newton analysis in Henlev et al.l 
(|2007D . We attribute this discrepancy to previously un- 




FlG. 1. — The shadowing filament used for our observations, 
shown in Galactic coordi nates. Grayscale: R OSAT All-Sky Sur- 
vey R1-I-R2 intensity ( S nowden et al.l liggTI) . Contours: IRAS 
100-micron intensity tS chleeel et al.lll998l ). Yellow squares: Our 
Suzaku pointing directions. 

recognized SWCX contamination in the XMM-Newton 
spectra, which means that SWCX contamination can oc- 
cur at times when the solar wind flux measured by A CE 
is low and does not show flares. 

This paper is organized as follows. The Suzaku obser- 
vations and data reduction are described in ^ The anal- 
ysis of the spectra using multicomponent spectral models 
is described in ^ The discrepancy between the Suzaku 
results and the XMM-Newton results is discussed in 2) 
This discrepancy is due to the presence of an additional 
emission component in the XMM-Newton spectra, which 
we also describe in [jH In f|5]we measure the total intensi- 
ties of the oxygen lines in our Suzaku and XMM-Newton 
spectra. We concentrate on these lines because they are 
the brightest in our spe ctra, and are a major com ponent 
of the 3/4-keV SXRB (jMcCammon et alJ [20021 ) . In g6] 
we present a simple model for estimating the intensity of 
the oxygen lines due to SWCX, which we compare with 
our observations. We discuss our results in SjTl and con- 
clude with a summary in SjSj Throughout this paper we 
quote la errors. 

2. OBSERVATIONS AND DATA REDUCTION 

Both of our observations were carried out in early 2006 
March. The details of the observations are shown in Ta- 
ble dJ In the following, we just analyze data from the 
back-illuminated XISl chip, as it is more sensitive at 
lower energies than the three front-illuminated chips. 

Our data were initially processed at NASA Goddard 
Space Flight Center (GSFC) using processing version 
1.2.2.3. We have carried out further processing and fil- 
tering, using HEAsoft^ v6.1.2 and CIAO^ v3.4. We first 
combined the data taken in the 3x3 and 5x5 observa- 
tion modes. We then selected events with grades 0, 2, 
3, 4, and 6, and cleaned the data using the standard 
data selection criteria given in the Suzaku Data Reduc- 
tion Guide'^. We excluded the times that Suzaku passed 
through the South Atlantic Anomaly (SAA), and also 
times up to 436 s after passage through the SAA. We 
also excluded times when Suzaku^ s line of sight was el- 
evated less than 10° above the Earth's limb and/or was 

^ http:/ /hcasarc. gsfc.nasa.gov/lhcasoft 
^ http://cxc.harvard.edu/ciao 

http: / /suzaku. gsfc.nasa.gov/docs/suzaku/analysis/abc/abc. html 
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TABLE 1 

Details of Our Suzaku Observations 





Observation 
ID 


/ 

(deg) 


b 

(deg) 


Start time 
(UT) 


End time 
(UT) 


Usable exposure 
(ks) 


Off filament 
On filament 


501001010 
501002010 


278.71 
278.65 


-47.07 
-45.30 


2006-03-01 16:56:01 
2006-03-03 20:52:00 


2006-03-02 22:29:14 
2006-03-06 08:01:19 


55.6 
69.0 



less than 20° from the bright-Earth terminator. FinaUy, 
we excluded times when the cut-off rigidity (COR) was 
less than 8 GV. This is a stricter criterion than that in 
the Data Reduction Guide, which recommends exclud- 
ing times with COR < 6 GV. However, the higher COR 
threshold helps reduce the particle background, and for 
observations of the SXRB one desires as low a particle 
background as possible. The COR threshold that we 
use has been used for other Suzaku observatio ns of the 
SXRB (|Fuiimoto et al.l[2007l : ISmith et al.ll20f)l . Finally, 
we binned the 2.5-8.5 keV data into 256-s time bins, and 
used the CIAO script analyze_ltcrv . si to remove times 
whose count-rates differ from the mean by more than 3(t. 
The resulting cleaned XISl images are shown in Figure^ 

We found that automated point source detection soft- 
ware did not work well on our Suzaku images, possibly 
because of Suzaku^ s rather broad point sprea d function 
(half-power diameter ~ 2'; iMitsuda e t al. '20071. There- 
fore, to remove sources which could contaminate our 
SXRB spectra, we used data from the pre-release of the 
second XMM-Newton Serendipitous Source Catalogue^ 
(2XMMp) to determine the locations of sources in our 
Suzaku fields of view. We excluded sources in the cat- 
alog with 0.2-4.5 keV fluxes above 5 x 10"^'' erg cm~^ 
s~^. For sources with fluxes in the range 5-30 x lO"^** 
erg cm~^ we excluded regions of the detector that 
were within 1.5' of the source. For brighter sources we 
excluded regions within 2.5' of the source. The excluded 
regions are shown in Figure [2l 

We extracted spectra from the full XISl field of view, 
excluding the above-mentioned sources. This field of 
view includes the parts of the chip illuminated by the 
Fe-55 calibration source. However, this will not affect 
our analysis, as we are only interested in X-ray ener- 
gies below those of the lines produced by the calibra- 
tion source. In our analysis we ignore channels 500-504 
(1.825-1.84325 keV), as there is artificial structure in 
the spectra at these energies^, due to a non-continuous 
conversion from pulse height amplitude (PHA) values to 
pulse invariant (PI) values near the Si K edge^. We 
binned up our spectra so that there were at least 25 
counts per bin. 

We constructed a particle background spectrum for 
each observation using tools available from the Suzaku 
Guest Observer Facility at GSFC^. These tools generate 
a background spectrum from night-Earth observations 
with the same distribution of cut-off rigidities as our 
observations. The background spectra were subtracted 
from the corresponding source spectra. Within the en- 
ergy range that we analyzed, there are instrumental lines 
in the background spectra at 1.49, 1.74, and 2.12 keV, 

^ http:/ /xminssc- www. star. le.ac.uk/newpages/xcat _public_2xmmp.html 
^ http:/ /suzaku. gsfc.nasa.gov/docs/suzaku/processing/ vl223.html 
^ http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/about/ug2/dotani.pdf 
http: / /suzaku. gsfc.nasa.gov/docs/suzaku/analysis/xisbgdO. html 



due to Al, Si, and Au, respectively. We found that 
these lines were not always accurately subtracted from 
our source spectra, leading to spurious features in the 
background-subtracted spectra. We therefore removed 
these three lines from the background spectra, and inter- 
polated the surrounding background continuum across 
this gap. In order to model these lines, we included three 
extra Gaussians in the models we fitted to our spectra. 

We generated redistribution matrix files (RMFs) using 
the tool mkxisrmf , and ancill ary response files ( ARFs) 
using the tool xissimarf gen ([Ishisaki et al.|[2007l ). This 
latter tool takes into account the spatially varying con- 
tamination on the optical blocking filters of the XIS sen- 
sors, which reduces th e detector efficiency at low energies 
(jKovama et al.l [20071 ). For the ARF calculation we as- 
sumed a uniform source of radius 20'. When generating 
the spectral response files, we used the set of calibration 
database (CALDB) files released on 2007 Jan 31. 

3. SUZAKU SPECTRAL ANALYSIS 
3.1. Description of Model 

We fit a model to our spectra which consists of compo- 
nents corresponding to the LB, the Galactic halo, and the 
extragalactic background due to unresolved active galac- 
tic nuclei. For the LB we used a single thermal plasma 
model in coUisional ionization equilibrium (CIE), while 
for the Galactic halo we used two equilibrium thermal 
plasma components. The use of a two-temperature halo 
model is well e stablished from analyses of the ROSAT 
All-S ky Survey (|Kuntz fc Snowdenll200(il : ISnowden etlHI 
|2000() . and is required to get a good fit to our data 
(see §3.2|, below). For the extragalactic background we 
used a power-law whose photon in dex was frozen at 1.46 
(jChen. Fabian k. Gendreaul [l997l ) . The LB component 
was not subject to any absorption, whereas the remain- 
ing components were. The temperatures of the thermal 
emission components and the normalizations of all four 
emission components were free parameters. 

We fit the above-described model simultaneously to 
our on- and off-filament 0.3-5.5 keV Suzaku spectra. 
Except for the adopted absorbing column densities, all 
of the model parameters were assumed to be the same 
for both spectra. The difference in absorbing column 
helps us to separate out the foreground emission (due to 
the LB) from the background emission (d ue to the halo 
aiid_the extragalactic background) . As in iHenlev et al.l 
(2007), we obtained the absorbing column densities for 
our observation directions from the IRAS 100-micron 
intensities for these directions (7.10 [on] and 1.22 [off] 
MJy sr-i; ISchlegel. Finkbeiner. fc DavisI Il998l ). These 
were converted to column den si ties u sing the conver- 
sion relation in ISnowden et al.l (|200Cl( ). The resulting 
on- and off-filament column densities are 9.6 x 10^" and 
1.9 X 10^" cni^^, respectively. 

We also include data from the ROSAT All-Sky Survey 
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Fig. 2. — Cleaned and smoothed Suzaku XISl images in the 0.3-5 keV band for our on-filament {lejt) and ofT-filament (right) observations. 
The data have been binned up by a factor of 4 in the detector's x and y directions, and then smoothed with a Gaussian wrhose standard 
deviation is equal to 1.5 times the binned pixel size. The particle background has not been subtracted from the data. The red circles 
outline the regions that were excluded from the analysis (see text for details). 



in our fit (|Snowden et al.lll997D . The Rl and R2 count- 
rates help constrain the model al lower energies (below 
~0.3 keV), while the higher channels overlap in energy 
with our Suzaku spectra. We extracted the ROSAT spec- 
tra from 0.5° radius circles centered on our two Suzaku 
pointing directions using the HEASARC X-ray Back- 
ground Tool^ v2.3. 

The s pectral analysi s was carried out using XSPEC^ 
vll.3.2 (|Arnaudl |1996( ). For the thermal plasma com- 
ponents, we used th e Astrophysical P lasma Emission 
Code (APEC) vl.3.1 (jSmith et aLllMll ) for the Suzaku 
data and the ROSAT R4-7 bands, a nd the Raymond- 
Smith code ([Raymond fc SmithI 119771 and updates) for 
the ROSAT Rl-3 bands. For a given model component 
(i.e., the LB or one of the two halo components), the tem- 
perature and normalization of the ROSAT Raymond- 
Smith model are tied to those of the corresponding 
Suzaku APEC model. We chose to use the Raymond- 
Smith code for the lower-energy ROSAT channels be- 
cause APEC's spectral calculations below 0.25 keV are 
inaccurate, due to a lack of data on transitions from L- 
sheh ions of Ne, Mg, Al, Si, S, Ar, and Ca^". As the 
upper-hmit of the ROSAT Rl and R2 bands is 0.284 keV, 
and the R3 band also includes such low-energy photons 
(jSnowden et al.lll997D . APEC is not ideal for fitting to 
these energy bands. 

For the absorption, we used the XSPEC 
phabs model, which uses cross-sections from 
iBaluciiiska-Church fc McCammonI |)1992D . 
with an updated He cross-section from 



Yan. Sa deghpour. fc Dalgarnol 
Henlev et al." (2003) 



1 



Following 
interstellar chemical 



^ , we used the 

abundance table from iWilms. Allen, fc McCravl (|2000D . 
For many astrophysically abundant elements, these 
abundances are lo wer than those in the widely used solar 
abundance table of lAnders fc Gre vesse (1989') . However, 
recently several elements' solar photosphc ric abundance s 
have been revised downwards (Asplu nd et al.l [20051) . 
and are in good agreement with the iWilms et al.l ()2000D 



http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/xraybg/xray bg.pl 
http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/xspecll 
^ http: / /cxc. harvard.edu/atomdb/issucs_cavcats. html 



abu ndances. The r efore, like Henl ev et al.l (|2007f ) we take 
the IWilms et al.l ()2000[ ) interstellar abundances to be 
synonymous with solar abundances. 

As noted in !j2l we included three Gaussians to model 
the Suzaku instrumental lines from Al, Si, and Au. The 
parameters of these lines were independent for the two 
Suzaku observations. 

During the course of the spectral modeling, it became 
apparent that there is a discrepancy in the normalization 
of the extragalactic power-law (EPL) between the two 
Suzaku spectra. For a photon index of 1.46, the EPL 
normalizations for the 2.0-5.5 keV energy range are ^^11 
(on-filament) and ~8 (off) photons cm~^ s~^ sr~^ keV~^ 
at 1 keV (cf. the ex pected value is ■- -^lO photons cm~^ 



expe 

s^i sr^i keV"^ e.g.. lChen et al.lll997[ ). We therefore al- 
lowed the normalization of the EPL to differ for the on- 
and off-filament spectra; however, the normalizations of 
the thermal plasma components were still constrained to 
be equal for the on- and off-filament spectra. We experi- 
mented with other methods for dealing with this discrep- 
ancy, but none of the results were significantly affected. 
We discuss this discrepancy further in i i7.3l 

3.2. Results 

Our on- and off-filament Suzaku spectra are shown 
in Figure [3l along with the best-fitting multicomponent 
spectral model described in the previous section. The 
model parameters are presented in Table [7| (Model 1). 
Overall, the model gives a good fit to the data: = 
734.24 for 703 degrees of freedom. However, the fit to 
some of the ROSAT bands is rather poor, as shown in 
Figure [H This could be due to a discrepancy in the 
effective area calibration between ROSAT and Suzaku. 
This discrepancy may be due to the uncertainty in the 
amount of contaminant on the XISl optical blocking fil- 
ter. For example, if the true amount of contaminant 
is larger than the amount given by the contamination 
database (CALDB) contamination model, the calculated 
effective area will be larger than the true effective area, 
and the resulting model normalizations will be too small. 

We investigated this possibility by adding a vphabs 
absorption component to our model, to model any con- 
tamination over and above that already included in the 



OBSERVING CHARGE EXCHANGE WITH SUZAKU AND XMM-NEWTON 



5 



On filament 



Local Bubble 

Halo (cold) 
Halo (hot) 
Extragalactic 

lental lines 

Total 



0.01 



0.001 




Off filament 



Local Bubble 

Halo (cold) 
Halo (hot) 
Extragalactic 

lental lines 

Total 



0.01 



0.001 



0.3 



1 

channel energy (keV) 




channel energy (keV) 



1000 



100 



Fig. 3. — Our observed on-filament {left) and off-filament {right) Suzaku spectra, with the best-fitting model obtained by fitting jointly 
to the Suzaku and ROSAT data (Model 1 in Table [7]l. The gap in the Si K instrumental line is where channels 500—504 have been removed 
from the data (see [[2]l. 

tion given by the CALDB contamination model, which 
is Nc = 3.1 X 10^^ cm"2 the center of the XISl chip 
(from the CALDB file ae_x;il_contami_20061016 .fits). 
This correction does not seem unreasonable, given that 
the systematic uncertainty on the contaminant thickness 
is -0.5 X lO^s cm^2 " 

We investigated whether or not a LB + two- 
temperature halo model is justified by fitting a model 
without a LB component to our data, and by fitting a 
LB + one-temperature halo model to our data. We also 
included a power-law to model the extragalactic back- 
ground. Both these models gave poor fits to the data: 
= 1068.63 and 870.51, respectively, for 704 degrees of 
freedom (Models 3 and 4 in Table [7]). In both cases, the 
large values of are mainly due to very poor fits to the 
ROSAT Rl and R2 data. 

Because CIE models give good fits to the data, we did 
not investigate non-equilibrium ionization models. How- 
ever, it should be noted that there are a number of indi- 
vidual features which are not well fit. In the off-filament 
spectrum, there is excess emission at ~0.5 keV. This 
feature may be too narrow to be an emission line, but 
if it is it is most likely due to N VII Lya. The Ne ix 
Hea emission at ~0.9 keV is poorly fit in both spec- 
tra - in both spectra the observed emission is in excess 
of the model. Finally, there is excess emission in the 
on-filament spectrum at ~1.3 keV, just below the Al K 
instrumental line at 1.49 keV. This is most likely due to 
Mg XI Hea. This emission may also be present in the 
off-filament spectrum, although it does not show up in 
the residuals. Note that the Gaussian representing the 
Al K line in the off-filament spectrum is broader and at 
a lower energy than that in the on-filament spectrum. It 
is possible that in the off-filament spectrum the data are 
not good enough to distinguish clearly the Mg xi Hea 
line and the instrumental Al K line, and that the Gaus- 
sian is in fact fitting to both lines, which would both 
broaden it and lower its energy. 
Table [7] also shows some variants of the above model. 




R3 R4 R5 
ROSAT band 

Fig. 4. — The on-filament {dashed) and off-filament {solid) 
ROSAT All-Sky Survey spectra, compared with Model 1 from Ta- 
ble [7] For clarity the individual model components have not been 
plotted. 



CALDB. This model component attenuated the emis- 
sion from the LB, halo, and extragalactic background 
for the Suzaku spectra only. We adjuste d the model 
oxyg en abundance to give C/0 = 6 (Ko vama et al.l 
|2007() . and set the abundances of all other elements to 
zero. The results of this model are given as Model 2 
in Table [T] Figures [5] and [6] show this model com- 
pared with the Suzaku and ROSAT spectra, respec- 
tively. One can see from Figure [6] that the fit to the 
ROSAT data is greatly improved. The model implies a 
column density of carbon atoms, Nc = (0.28 ± 0.04) x 
10^^ cm^'^, in addition to the amount of contamina- 



http: / /www. astro, isas.ac.jp/suzaku/proccss/cavcats/cavcats_xrtxis08. html 



6 



HENLEY AND SHELTON 



On filament Local Bubble 



Halo (hot) 
Extragalactic 

nental lines 

Total 



0.01 



0.001 




0.01 



0.001 




1 


1 ■ 1 _ 







0.3 



1 

channel energy (keV) 



0.3 



1 

channel energy (keV) 



Fig. 



As Figure [3l but with a vphabs component included for the Suzaku spectra to model XIS contamination above that included 



in the CALDB (Model 2 in Table [7] see f|32]for details). 
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Fig. 6.— As Figure m but for Model 2 from Table [7] 



In Model 5 we fit exactly the same model to the Suzaku 
data alone. Without the ROSAT data, we cannot con- 
strain the LB temperature Tlb- We therefore fix it at 
the value determined in Model 2: Tlb = 10^-^^ K. We 
also fix Nq for the vphabs contamination component 
at the Model 2 value: iVc = 0.28 x lO^^ cm'^. The 
best-fitting model parameters are in very good agree- 
ment with those obtained by fitting jointly to the Suzaku 
and ROSAT spectra (compare Models 2 and 5). We can 
only get consistent results between the Suzaku- ROSAT 
joint fit and the fit to just the Suzaku data by using a 
two-temperature model for the halo. However, note that 
the Model 5 LB emission measure is consistent (within 
its errorbar) with zero. This is not to say that our data 
imply that there is no LB at all: as noted above, we need 
a LB component and two halo components to get a good 
joint fit to the Suzaku and ROSAT data. Instead, the 
Model 5 results imply that our Suzaku data are consis- 



tent with the LB not producing significant emission in 
the Suzaku band (i.e., E > 0.3 keV). 

Also shown in Table [7] are the results of fitting our 
model (without any LB component) to the on- and off- 
filament Suzaku spectra individually (Models 6 and 7, 
respectively). As has already been noted, the normaliza- 
tion of the extragalactic background differs significantly 
between the two spectra. However, the plasma model pa- 
rameters for the individual spectra are in good agreement 
with each other. These results seem to justify allowing 
the extragalactic normalization to differ between the two 
spectra while keeping all other model parameters equal 
for the two spectra. 

3.3. Chemical Abundances in the Halo 

We investigated the chemical abundances of the X-ray- 
emissive halo gas by repeating the above-described mod- 
eling, but allowing the abundances of certain elements in 
the halo components to vary. In particular we wished to 
investigate whether or not varying the neon and magne- 
sium halo abundances improved the fits to the Ne ix and 
Mg XI features noted above. We also allowed the abun- 
dance of iron to vary, as iron is an important contributor 
of halo line emission to the Suzaku band. 

For this investigation we just fit to the Suzaku data, 
fixing the LB temperature at log(rLB/K) = 5.98, and 
fixing the carbon column density of the vphabs contam- 
ination component at A^c = 0.28 x 10^^ cm~^. As we 
could not accurately determine the level of the contin- 
uum due to hydrogen, we could not measure absolute 
abundances. Instead, we estimated the abundances rel- 
ative to oxygen by holding the oxygen abundance at its 
I Wilms et alj ()200C1| ) value, and allowing the abundances 
of neon, magnesium, and iron to vary. Both halo com- 
ponents were constrained to have the same abundances. 

The best-fitting temperatures and emission measures 
of the various model components are presented as 
Model 8 in Table [7l and the abundances are presented 
in Table [21 The best-fitting model parameters are not 
significantly affected by allowing certain elements' abun- 
dances to vary (compare Model 8 with Model 5). Iron 
does not seem to be enhanced or depleted relative to oxy- 
gen in the halo. Neon and magnesium both appear to be 
enhanced in the halo relative to oxygen, which is what 
one would expect from Figures [3] and [5l as the models 
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TABLE 2 
Halo Abundances 

Element Abundance 



O 
Fe 



1 (fixed) 
1.8 ±0.4 
^+3.5 

-2.8 



4.6"' 



1.2 



+0.4 

0.; 



Note. — Abun- 
dances are rela ti ve to 
the W ilms et al.l ll2000h 
interstellar abundances: 
Ne/O = 0.178, Mg/O = 
0.051, Fe/O = 0.055. 
^ These enhanced abun- 
dances may be an arte- 
fact of SWCX contami- 
nation; see i|7.2l 



shown in those figures underpredict the neon and mag- 
nesium emission. 

We discuss these results in §ij7?2land In i i7.2l we 

discuss the possibility that the enhanced neon and mag- 
nesium emission is in fact due to SWCX contamination 
of these lines, rather than being due to these elements 
being enhanced in the halo. On the other hand, in i j7.5l 
we discuss the implications of neon really being enhanced 
in the halo with respect to oxygen and iron. 

4. COMPARING THE SUZAKU AND XMM-NEWTON 
SPECTRA 

For comparison, Table[7]also contains the results of the 
analysis of the XMM- Newton spectra from the same ob- 
servation directions bv lHenlev et al.l ()2007f ). The Model 9 
results are taken directly from thei r "standard" model . 
However, it should be noted that iHenlev et al.l ()2007f ) 
used APEC to model all of their data, whereas in the 
analysis described above we used the Raymond-Smith 
code to model the ROSAT Rl-3 bands. We have 
therefore reanalyzed the XMM-Newton -\- ROSAT spec- 
tra, this time using the Raymond-Smith code to model 
the ROSAT Rl-3 bands, and using APEC to model 
the ROSAT R4-7 bands and the XMM-Newton spec- 
tra. This new analysis allows a fairer comparison of 
our XMM-Newton results with our Suzaku results. The 
XMM-Newton spectra we analyzed are identical to those 
used by Henley ct al. (2007) - see that paper for details 
of the data reduction. We added a broken power-law to 
the model to take into account soft-proton contamination 
in the XMM-Newton spectra. This broken power-law 
was not folded through XMM-Newton's effective area, 
and was allo wed to differ for t he on- and off- filament 
datasets (see iHenlev et al1l2007li ). The presence of this 
contamination means we cannot independently constrain 
the normalization of the extragalactic background. We 
therefore freeze the on- and off-filament normalizations 
at the S'Mzafcw-determined values. 

The results of this new analysis are presented as 
Model 10 in Table [71 As can be seen, there is poor 
agreement between the best-fit parameters of the Suzaku 
+ ROSAT model (Model 2) and the XMM-Newton + 
ROSAT model (Model 10). We believe this discrepancy 
is due to an extra emission component in the XMM- 
Newton spectra. In Figure [7| we plot the differences 
between the XMM-Newton spectra and our best-fitting 
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Fig. 7. — The excesses in our on-filament (top) and off-filament 
(bottom) XMM-Newton spectra over our best-fitting model to the 
Suzaku -\- ROSAT data (Tabled Model 2). The gap in the data 
between 1.4 and 1.9 keV is where two bright instrumental lines 
have been removed. 



Suzaku + ROSAT model (Model 2 from Table [T]). To 
our best-fitting Suzaku -\- ROSAT model we have added 
a broken power-law to model the soft-proton contami- 
nation in the XMM-Newton spectra. The parameters 
of this broken power-law are frozen at the values de- 
termined from the fitting to the XMM-Newton spectra 
described in the previous paragraph. The on-filament 
XMM-Newton spectra show excess line emission at ~0.57 
and ~0.65 keV, most likely due to O vii and O viii, re- 
spectively. The features in the off-filament spectra are 
not as clear. However, there appears to be excess O vii 
emission in the MOSl spectrum, and excess emission at 
~0.7 keV (of uncertain origin) and '^0.9 keV (proba- 
bly Ne ix) in the M0S2 spectrum. We can estimate 
the significance of the excess emission by calculating 
for the XMM-Newton data compared with the Suzaku 
+ ROSAT model. We concentrate on the excess oxy- 
gen emission and calculate for the 0.5-0.7 keV energy 
range. We find — 106.28 for 24 degrees of freedom for 
the on-filament spectra, and = 43.03 for 22 degrees of 
freedom for the off-filament spectra. These correspond 
to probabihties of 2.5 x 10~^^ and 0.0047, respectively, 
implying that the excesses are significant in both sets of 
spectra at the 1% level. 

We measure the intensities of the extra oxygen emis- 
sion by fitting (5-functions at E — 0.570 keV and 
0.654 keV to the excess spectra in Figure [T] We fit 
these (5-functions simultaneously to the on- and off- 
filament XMM-Newton excess spectra. The intensities 
of the excess oxygen emission in the XMM-Newton spec- 
tra over the best-fitting Suzaku + ROSAT model are 
3.8 ± 0.5 L.U. (O vii) and 1.4 ± 0.3 L.U. (O viii). 

We believe that this excess oxygen emission is due to 
SWCX contamination in our XMM-Newton spectra. As 
noted in the Introduction, this was not taken into ac- 
count in the original analysis of the XMM-Newton spec- 
tra. This is because the solar wind flux was steady and 
slightly below a verage during t he XM M-Newton obser- 
vations, leading iHenley et al.l (|2007t ) to conclude that 
SWCX contamination was unlikely to be significant. We 
discuss the SWCX contamination in our spectra further 
in P 

5. MEASURING THE OXYGEN LINES 
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As well as using the above-described method to sep- 
arate the LB emission from the halo emission, we mea- 
sured the total intensities of the O vii complex and O viii 
line at ^^0.57 and ~0.65 keV in each spectrum. These 
lines are a major component of the SXRB, accounting 
for the majority of the observed ROSAT R4 diffuse back- 
ground that is not due to resolv ed extragalactic discrete 
sources (jMcCammon et aLll2002f) , and are easily the most 
prominent lines in our Suzaku spectra. 

To measure the oxygen line intensities, we used a model 
consisting of an absorbed power-law, an absorbed APEC 
model whose oxygen abundance is set to zero, and two 
(5-functions to model the oxygen lines. As in the pre- 
vious section, the power-law models the extragalactic 
background, and its photon index was frozen at 1.46 
(|Chen et al.lll997f) . The APEC model, meanwhile, mod- 
els the line emission from elements other than oxygen, 
and the thermal continuum emission. The absorbing 
columns used were the same as those used in the ear- 
lier Suzaku analysis. As with our earlier analysis, we 
multiplied the whole model by a vphabs component to 
model the contamination on the optical blocking filter 
which is in addition to that included in the CALDB con- 
tamination model (see t j3.2p . We fix the carbon column 
density of this component at 0.28 x 10^^ cm~^ (Table[7l 
Model 2). We fit this model simultaneously to our on- 
and off-filament spectra. However, all the model param- 
eters were independent for the two directions, except for 
the oxygen line energies - these were free parameters 
in the fit, but were constrained to be the same in the 
on- and off-filament spectra. We used essentially the 
same method to measure the oxygen line intensities in 
our XMM-Newton spectra. However, we did not use a 
vphabs contamination model, and, as before, we added 
a broken power-law to model the soft-proton contami- 
nation. Table [3] gives the energies and total observed 
intensities of the O vii and O viii emission measured 
from our Suzaku and XMM-Newton spectra. 

We can use the difference in the absorbing column for 
the on- and off-filament directions to decompose the ob- 
served line intensities into foreground (LB + SWCX) and 
background (halo) intensities. If Ifg and /halo are the in- 
trinsic foreground and halo line intensities, respectively, 
then the observed on-filament intensity /on is given by 

/on = /fg + e~'^°°/halo, (1) 

where Ton is the on-filament optical depth at the energy 
of the line. There is a similar expression for the observed 
off-filament intensity log, involving the off-filament op- 
tical depth Toff. These expressions can be rearranged to 
give 

e ""/on ^ °^^/off /(^\ 
T ^on /off /„\ 

For the purposes of this decomp osition, we use the 
iBalucihska-Church fc McCammonI (|T992i') cross-se ctions 
(wi th an updated He cr oss-section; [Yan et al]ll998f ) with 
the lWilms et al.l ()2000D interstellar abundances. We use 
the cross-sections at the measured energies of the lines. 
For the Suzaku measurements, the cross-sections we use 
are 7.17 x 10"22 cm^ for O vii {E = 0.564 keV) and 
4.66 X 10-22 for O viii {E = 0.658 keV). For the 



XMM-Newton O vii emission we use a cross-section of 
7.03 X 10^22 ^ 0.568 keV). We cannot decom- 

pose the XMM-Newton O viii emission because the on- 
filament O VIII line is brighter than the off-filament line. 
This gives rise to a negative halo intensity, which is un- 
physical. 

The results of this decomposition are presented in Ta- 
ble m Note that the foreground oxygen intensities mea- 
sured from the Suzaku spectra are consistent with zero. 
This is consistent with our earlier finding that the Suzaku 
spectra are consistent with there being no local emis- 
sion in the Suzaku band (see §3.2p . The difference be- 
tween the foreground O Vii intensities measured from 
our XMM-Newton and Suzaku spectra is 5.1 ± 3.1 L.U.. 
This is consistent with the O vii intensity measured from 
the excess XMM-Newton emission over the best-fitting 
Suzaku + ROSAT model (3.8 ± 0.5 L.U.; see ©. The 
halo O VII intensities measured from our XMM-Newton 
and Suzaku spectra are consistent with each other. This 
is as expected, as we would not expect the halo intensity 
to significantly change in '^4 years. 

6. MODELING THE SOLAR WIND CHARGE EXCHANGE 
EMISSION 

In §Sj4] and [5] we presented evidence that our XMM- 
Newton spectra contain an extra emission component, in 
addition to the components needed to explain the Suzaku 
spectra. In particular, the O Vii and O viii emission are 
enhanced in the XMM-Newton spectra. We attribute 
this extra component to SWCX emission, as it seems 
unlikely to be due to a change in the Local Bubble or halo 
emission. This extra component helps explain why our 
XMM-Newton and Suzaku analyses give such different 
results in Table [T] 

Previous observations of SWCX have found that in- 
creases in the SWCX emission are associated with en- 
hancements in the solar wind, as measured by ACE. 
These enhancements consist of an increase in the pro- 
ton flux, and may also include a shi ft in the ionization 
balance to higher ion ization stages ( Snowden et al.l[2004t 
iFunmoto et al.ll2007h . In §6.1) we present a simple model 
for heliospheric and geocoronal SWCX emission, and use 
contemporaneous solar wind data from the ACE and 
WIND satellites to determine whether or not the ob- 
served enhancement of the oxygen lines in the XMM- 
Newton spectra is due to differences in the solar wind 
between our two sets of observations. 

In addition to the variability associated with solar wind 
enhancements, the heliospheric SWCX intensity is also 
expected to vary during the solar cycle, due to the dif- 
ferent states of the solar wind at solar maximum and 
solar minimum (jKoutroumpa et aT]|2006D . As our two 
sets of observations were taken ~4 years apart, at differ- 
ent points in the solar cycle, in §6. 21 we examine whether 
the SWCX intensity variation during the solar cycle can 
explain our observations. 

6.L A Simple Model for Heliospheric and Geocoronal 
SWCX Emission 

6.1.1. The Basics 

A SWCX line from a X+" ion of element X results 
from a charge exchange interaction between a X+("+^^ 
ion in the solar wind and a neutral atom. The intensity 
of that line therefore depends on the density of X+("+^^ 
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TABLE 3 

Observed Total Oxygen Line Intensities 



Satellite 



Dataset 



Energy 
(keV) 



O vn 

intensity 
(L.U.) 



Energy 
(keV) 



O VIII 

Intensity 



(L.U.) 



Suzaku 



XMM-Newton 



On filament 
Off filament 
On filament 
Off filament 



0.564 6.51 
10.53 

0.568 10.65 
13 



+0.37 
0.45 
+0.68 
0.55 
+0.77 
0.82 
+ 1.58 
-1.49 



0.658 



0.656 



2.54 
3.21 
3.91 
2.81 



+0.22 
0.28 
+0.25 
0.36 
+0.32 
0.20 
+0.58 
0.60 



Note. — L.U. = photons cm ^ s -"^ sr ^. 



TABLE 4 

Foreground and Halo Oxygen Line Intensities 





O VII 


O VIII 


Satellite 


(L.U.) 


Jhulo 

(L.U.) 


Jfg ^halo 

(L.U.) (L.U.) 


Suzaku 

XMM-Newton 


-'-■-'--1.4 
f\ 9+2.8 
"■^-2.9 


°-°-4.6 


1.0 ±1.1 2At\\i 


Note. — L.U. 


= photons 


cm~^ s~ 


■1 sr-i. 



ions in the solar wind, nx+cn+i), and on the density of 
ne utral atoms, Un. The line intensity, J, can be written 
as (|Cravensll200"0l: IWargelin et al.ll2004D 



/ = [ cryn^ j "•x+'"+i) j f !!2i \ n^vjUswdl, (4) 
in J \ nx J yniij 

where a is the cross-section for a charge exchange reac- 
tion between a X+("+^) ion and a neutral, y is the yield 
of the particular line of interest, rix+cn+i) /nx is the ion 
fraction of the X+("+^) ion, nx/nn is the abundance of 
element X in the solar wind, and rig^ and the 
proton density and speed of the solar wind. If, for ex- 
ample, we wish to calculate the intensity of the O vii 
emission at ^^0.57 keV, the "line" of interest is the blend 
of n = 2 ^ 1 transitions of 0+^ ions, and the relevant 
ion fraction to insert in the integrand in equation ^ is 
that of 0+^ ions. As was stated in the Introduction, 
there are contributions to the SWCX emission from so- 
lar wind ions interacting with neutral interstellar H and 
He atoms distributed throughout the solar system (hclio- 
spheric emission) , and with neutral H atoms in the outer 
reaches of the Earth's atmosphere (geocoronal emission). 
For the heliospheric emission, one must calculate the con- 
tributions due to interactions with H and He separately, 
and add them. This is because the cross-sections and 
yields for these interactions are different, as are the den- 
sities of H and He in the heliosphere. 

In the following, we estimate the intensities of the O vii 
and O VIII emission due to heliospheric and geocoronal 
SWCX emission for our XMM-Newton and Suzaku ob- 
servations. For this purpose, we adopt simple models for 
the density of the neutral atoms, and use contemporane- 
ous solar wind data from the ACE and WIND satellites 
to insert in equation To calculate the O vii inten- 
sity, we use a = 3.40 x 10"^^ and 1.80 x 10"^^ cm^ 
for O"*"^ interacti ng via charge exchange with H and 
He, respectively (jKoutroumpa et aT] |2006[ ). The yield 
of O VII Ka emission from interactions between 0+^ 



and He is y = 0.86, where y is the sum of the contribu- 
tions of the resonance, intercombination, and fo rbidden 
hues (jKrasnopolskv. Greenwood, fc Stancim2004l ). Yield 
information is not available for 0+^ interacting with 
H. However, by definition y < 1, so we can calcu- 
late an upper-limit for the expected SWCX emission 
by setting y = 1. To calculate the O viii intensity, 
we use a = 5.65 x lO'^^ and 2.80 x 10'^^ cm^ for 
charge exc hange interactions of O^^ with H and He, re- 
spectively (|Koutroumpa et aT]|2006l ). The yield of the 
O VIII Lyg line due to ch arge exchange with He is 0.65 
(jKrasnopolskv et al.ll2004) . Again, yield information is 
not available for 0+^ interacting with H, so we calculate 
an upper limit on the expected intensity by setting y = 1. 

6.1.2. Solar Wind Data 

We use solar wind data from the ACE and WIND 
satellites, downloaded from the ACE Science Center^^ 
and CDAWeb^^, respectively. The solar wind data from 
around the times of our XMM-Newton and Suzaku obser- 
vations are shown in Figure [H Note the different ranges 
on the time axes - this is because the Suzaku point- 
ings were much longer than the XMM-Newton point- 
ings. Where possible, we use the proton density and 
wind speed (and hence proton flux) from the ACE Solar 
Wind Electron, Proton, and Alpha Monitor (SWEPAM). 
In Figure m we just plot good data (bad SWEPAM data 
are denoted by a value of —9999.9 in the density or speed 
column, and so cannot be plotted in Fig. [5]). As can 
be seen in Figures [DJf-h), there is a gap in the ACE 
SWEPAM dataset (plotted in black), which is almost 
exactly coincident with our on-filament Suzaku observa- 
tion. We use data from the WIND Solar Wind Experi- 
ment (SWE) to fill in this gap; these data are plotted in 
red in Figures [Hf-h) . In the WIND dataset, good data 
are denoted by a quality flag of 0. Unfortunately, for 
the times covered by Figures [8l^f-h), the WIND dataset 
contains no good data, and hence should not be used. 
However, unlike the ACE SWEPAM data, densities and 
speeds are quoted for bad data points, so the data can 
still be plotted. As can be seen in Figures E^f-h), at 
times when there is both WIND data and ACE data, 
the two satellites are in good agreement regarding the 
shape of the solar wind variations. This suggests that 
the WIND data should give a reasonably trustworthy 
picture of the solar wind proton flux in the gap in the 
ACE data. Figure [5] also shows oxygen ion data from the 

12 http://www.srLcaltcch.odu/ACE/ASC/lGvol2 
^'^ http://cdawcb.gsfc.nasa.gov 
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TABLE 5 

Solar Wind Parameters Used In SWCX Model 





nswo^ 


1/ b 






Satellite 


(cm"'') 


(km s~^) 






Suzaku 


8.3 


360 


0.13 


5.6 X lO""* 


XMM-Newton 


4.2 


420 


0.15 


5.6 X 10"* 



Solar wind proton density (from ACE and WIND). 
^ Solar wind speed (from ACE and WIND). 
Solar wind ion fraction (from ACE). 

So lar wind oxygen abundance (from ISchwadron fc CravensI 
[2000D 

ACE Solar Wind Ion Composition Spectrometer and So- 
lar Wind Ion Mass Spectrometer (SWICS/SWIMS). The 
data shown are two-hour averages. Here we only show 
good data: those with a quality flag of 0. 

The ACE data from the time of the XMM-Newton ob- 
servations show that the solar wind was very steady at 
that time, and that the proton flux and O"*"^ ion frac- 
tion were slightly below their typ ical values. As alrea dy 
mentioned, these observations led lHenlev et al.l ()2007D to 
conclude that SWCX contamination was likely to be low 
in the XMM-Newton spectra. However, the results of 
the previous sections imply that this is not the case, and 
that SWCX emission makes up most of the foreground 
oxygen emission in the XMM-Newton spectra. 

Around the time of the Suzaku observations, the solar 
wind was much more variable. In particular, the pro- 
ton flux rises by a factor of ~4 during the last third of 
the on-filament Suzaku observation. However, we find no 
increase in the soft X-ray (0.3-2.0 keV) count-rate asso- 
ciated with this increase in the proton flux. Furthermore, 
the results of our Suzaku spectral fitting do not signif- 
icantly change if we use an on-filament spectrum that 
excludes the final third of the observation, instead of a 
spectrum from the whole of the on-filament observations. 
These results suggest that there are not changes in the 
SWCX emission directly associated with the changes in 
the solar wind during the Suzaku observations. 

From the data shown in Figure [3 we obtain values of 
no+r/no, Usw, and Usw to insert into equation (jl]), in 
order to estimate the expected O Vli intensity due to 
SWCX during our XMM-Newton and Suzaku observa- 
tions. These values are shown in Tabled The ACE 
data we obtained do not include the 0+^ ion fraction, 
which is needed to estimate the expected O viii inten- 
sity. Instead, we use the ra tio 0+^/0+^ = 0.35 from 
ISchwadron fc Cravend ()2000D and our measured 0+'' ion 
fractions to estimate the O"*"^ ion fraction. In princi- 
ple we could also obtain no / from A CE. In prac- 
tice this would be measured from the He/0 ratio from 
SWICS/SWIMS. However, these data exhibit an unex- 
pected 6-month periodicity, meaning that they are un- 
reliable (K. D. Kuntz, private communication). Instead, 
we use a canonical sl ow solar wind oxygen abunda nce of 
no/nu = 5.6 x lO"" ()Schwadron fc Cra.vensl[20M) . 

6.1.3. Geocoronal Emission 

For the geocoronal emission, we assume that the neu- 
tral H density varies with geocentric distance Tg as 

25 cm~'^ is the 



exospheric neutra l H density at 10-Rei w here i?E is the ra- 
dius of the Earth (jWargehn et al.l[2"004[ ). We also assume 
that the solar wind ion density is zero inside th e magneto- 
spher e, and constant everywhere outside (War gelin et al.l 
|2004[ ). With these assumptions, equation (j4|) yields 



1 f°° 
-^o VII (geocoronal) = — Q!n„o»^swo"sw j 



/10i?E\' 



47r 



\ ' min / 



dra 



(5) 



where a — (7y{nQ+7 /no){no /nn) , Ugwo is the solar wind 
density at 1 AU, and r„iin is the distance from the Earth 
to the magnetopause. For simplicity, we have integrated 
equation (|4]) radially away from the Earth, rather than 
along the line of sight from the satellite. This approxima- 
tion is likely to have a bigger effect on the XMM-Newton 
result, as Suzaku is in a low-Earth orbit. In principle, a 
more accurate intensity could be obtained by integrating 
equation ((4]) numerically along the line of sight. How- 
ever, as the geocoronal emission is an order of magnitude 
fainter than the heliospheric emission (see below), the ap- 
proximation used in equation ([5]) is unlikely to affect our 
conclusions. 

For our observations, the angle between the line of sight 
and the direction towar d the Sun was approximately 80° . 
From Figure 7 in Warg ehn et all (2004) , we estimate that 
''min ~ 13i?E for this viewing angle. Given the other 
approximations in the model, a more accurate estimate 



6.1.4. Heliospheric Emission 

To estimate the heliospheric contribution to the SWCX 
emission, one needs to know the density of neutral H 
and He throughout the heliosphere. Neutral H and He 
are depleted in the inner solar system, due to photoion- 
ization and charge exchange (jCravend |2000( ). This ef- 
fect is stronger for H than for He. Although it is pos- 
sible to construct detailed models of the density distri- 
butions of interstellar ne utrals in the solar system (e.g., 
iKoutroumpa et all l20()6t see ^6.21 below), for this sim- 
ple estimate we follow ICravensI ()2000^ and use nn(r/i) = 
rino exp(— A/r/i), where rh is the heliocentric distance, 
rino is the interstellar density, and A is the attenuation 
length. We adopt — 0.15 cm~'^ and A = 5 AU 
for H, and Unp — 0.015 cm""^ and A = 1 AU for He 
(jCravensI IMifif? . We also assume that the solar wind 
is isotropic, so the solar wind density Ugw varies as 
nsw{rh) = nswoifo/fh)'^', as before, n^wo is the solar wind 
density a,t rh = tq = 1 AU. We integrate from the Earth 
to the edge of the heliosphere. As the distance to the he- 
liopause is much larger than A, we can replace the upper 
limit with infinity. If, for simplicity, we were to integrate 
radially away from the Sun, rather than along the line of 
sight from the Earth, equation ^ would yield 



^ovii(liehospheric) = -^-arino'^swo'^sw / 



Th 



47rA 



1 -e' 



-A/ro 



nr,{rg) = n„o(10i?E/?'g) , where n„ 



where a is defined as before. In practice, we integrate 
equation (|4|) numerically along the line of sight. This 
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Fig. 8. — Solar wind data from around the times of our XMM-Newton {left) and Suzaku {right) observations. Panels (a) and (f) show 
the solar wind proton flux, measured by the ACE SWEPAM experiment {black) and, for the Suzaku observations, by the WIND SWE 
experiment as well {red). The proton flux is the product of the proton density and the solar wind speed, shown in panels (b) and (c) for 
the XMM-Newton observations, and panels (g) and (h) for the Suzaku observations. Panels (d) and (i) show the O"'"'' ion fraction, and 
panels (e) and (j) show the ratio of to 0+^. The oxygen i on fraction data are fro m th e AC E SWICS/SWIMS experiment. The black 
dashed lines show typical values for the slow solar wind, from r Wargelin et al.l | |2004| ) and ISch wadron & Cravens (2000). The solar wind 
data have not been shifted to allow for the travel time from the spacecraft to the Earth {~4 ks). 



yields heliospheric intensities that are ^--^15% larger than 
those obtained from equation ([5]). 

We estimate the intensity of the O vii and O viii emis- 
sion due to geocoronal and heliospheric SWCX emission 
for our XMM-Newton and Suzaku observations using the 
yields and cross-sections from i i6.1.1l and the solar wind 
parameters from Table [5] The results are shown in Ta- 
ble [S 

6.2. The Heliospheric Emiss i on Mo del of 
\Koutroumva et ai\ i200di ) 



An alternatiy e ap proach, aeveiopea Dy 
iKoutroumpa et all ()2006[ ). separately considers the 
contributions of the slow and fast solar winds to the 
heliospheric SWCX emission, and how these contribu- 
tions change with direction and in relation to the solar 
cycle. They use a more sophisticated, self-consistent 
model for the distribution of interstellar neutrals within 
the solar system, taking into account the effects of 
gravity, radiation pressure, and losses due to charge 
exchange and photoionization. They also calculate 
self-consistently the distribution of ions in the solar 
wind, as this too will vary under the effects of charge 
exchange. Rather than using solar wind data from ACE 
for a particular time, as we did above, they adopt typical 
wind velocities, densities, ion fractio ns, and abundances. 

Mos t importantly for our purposes. IKoutroumpa et al.l 
(|2006f ) calculate heliospheric intensities at solar maxi- 



developed by 



mum and at solar minimum. This is relevant to our 
work because our XMM-Newton observations were taken 
in 2002 May, near the end of the last solar maxi- 
mum, whereas our Suzaku observations were taken in 
2006 March, at solar minimum. At the high ecliptic 
latitude of our observat ion directions (w —74°), the 
IKoutroumpa et all ()2006D model predicts higher O viii 
intensities at solar maximum than at solar minimum. 
This is due to a c hange in the state of the solar wind. 
At solar maximum, IKoutroumpa et al] (|2006( ) assume an 
isotropic "slow" solar wind, and at solar minimum they 
assume a slow solar wind in the solar equatorial zone, and 
a "fast" solar wind at heliographic latitudes above -1-20° 
and below —20° (the changes in the solar wind during 
the solar cycle are reviewed in lSmith et al.|[2003f) . The 
fast solar wi nd is in a lower ionization sta te than the slow 
solar wind (jSchwadron fc Cravens! [2000l ) . For example, 
the ion fractions of and 0+^ in the fast solar wind 
are 0% and 3%, against 7% and 20% for the slow so- 
lar wind. These di f ferent ion fractions explain why the 
IKoutroumpa et all ()2006D model gives O Viii intensities 
at high ecliptic latitudes that are lower at solar minimum 
than at solar maximum. At solar minimum a significant 
length of the line of sight passes through the fast solar 
wind, in which there is no O"*"^ and so no contribution to 
the O VIII emission. In contrast, at solar maximum there 
are contributions to the heliospheric O viii emission all 
along the line of sight to the heliopause. 
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TABLE 6 

Model and Observed SWCX Line Intensities 







Model SWCX intensity^* 


KoutroumDa et al. (2007) 


Observed SWCX 






Geocoronal 


Heliospheric 


TotaP 


model intensity 


intensity 


Satellite 


Ion 


(L.U.) 


(L.U.) 


(L.U.) 


(L.U.) 


(L.U.) 


Suzaku 


O VII 


0.28 


3.4 


3.7 


0.83 


< 3.3^^ 




O VIII 


0.16 


1.9 


2.0 


0.07 


< 3.2^^ 


XMM-Newton 


O VII 


0.19 


2.3 


2.5 


2.32 


3.8±0.5<* - 7.1 ±0.5" 




O VIII 


0.11 


1.3 


1.4 


0.92 


1.4±0.3<i -4.6 ±0.3" 



Note. — L.U. = photons cm s ^ sr \ 
^ Calculated using the model discussed in i|6.1l 

^ Calculated from the unrounded geocoronal and heliospheric intensities. 
2rT upper-limit on the foreground oxygen intensity from Table |4] 

Oxygen line intensity measured from the excess of the XMM-Newton emission over the best-fitting Suzaku 
ROSAT model (see gj. 
Sum of the XMM- Afewton-determined lower limit and the Suzaku-determined upper limit. 



In a follow-u p to iKoutroumpa et al] (|2006D , 
iKoutroumpa et al] (|2007f ) calculate heliospheric oxy- 
gen intensity for various XMM-Newton and Suzaku 
observations of the SXRB, including those discussed 
here. In Table [6] we give the O VII and O VIII inten- 
sities predicted by their "ground level" model for our 
observations. These "ground level" values are calcu- 
lated assuming the solar wind parameters are at their 
nominal values, without taking into account possible 
enhancements in the solar wind. 

6.3. Comparison with Observations 

As well as the SWCX intensities predicted by the mod- 
els discussed above, Table |6] contains the limits on the 
amount of SWCX oxygen emission observed by Suzaku 
and XMM-Newton. The upper limits determined from 
the Suzaku spectra are the 2a upper limits on the fore- 
ground oxygen intensities presented in Table ID We can- 
not make these upper limits any tighter, as we do not 
know a priori how much of the foreground emission is 
due to the LB, and how much is due to SWCX. In 2] we 
measured the intensities of the oxygen lines in the XMM- 
Newton spectra above the best- fitting Suzaku + ROSAT 
model. If we make the reasonable assumption that this 
excess emission is due to SWCX, the excess oxygen line 
intensities are lower limits on the SWCX oxygen emission 
observed by XMM-Newton. They are lower limits be- 
cause we cannot rule out the possibility that our Suzaku 
spectra (which we use as a baseline) contain some SWCX 
emission. We can obtain an upper limit on the SWCX 
oxygen emission observed by XMM-Newton by adding to 
this lower limit the upper limit on the amount of SWCX 
emission in the Suzaku spectra. 

There is poor agreement between the predictions of the 
model discussed in H6.1\ and the observed O VII line in- 
tensities. The model overpredicts the O vii intensity 
observed by Suzaku, and underpredicts that observed 
by XMM-Newton. As the discrepancies are in oppo- 
site senses for the two datasets, they cannot be reduced 
by altering input model parameters, such as the cross- 
sections, line yields, or the attenuation lengths A of neu- 
tral H and He in the heliosphere (see §6.1.4p . This is be- 
cause the same set of parameters is used to model both 
datasets, so reducing one discrepancy will increase the 
other. 



It is likely the model presented in §671] overpredicts 
the SWCX O VII intensity for Suzaku because we es- 
sentially assume an isotropic slow solar wind when in- 
tegrating equation (jH) to calculate the heliospheric in- 
tensity. This is not valid for observations taken at so- 
lar minimum (i.e., the Suzaku observations) as much 
of the line of sight passes through the fast solar wind, 
which produces less O vii and O viii emission because of 
the lo wer 0+^ and 0+ ^ ion fractions (|Koutroumpa et al.l 
l2006l ). However, the IKoutroumpa et aP (|2006l ) model 
does take into account the state of the solar wind at 
solar minimum, and one can see from Table [5] that 
it predicts lower heliospheric oxyge n intensities for the 
Suza ku observations than our model (|Koutroumpa et ahl 
l2007f) . The solar minimum oxygen in tensities predicted 
by the IKoutroumpa et al.l (|200'6l i2007f ) model are consis- 
tent with the Suzaku observations. 

As our model assumes an isotropic slow solar wind, 
the predicted intensities should be more applicable to 
observations taken at solar maximum (i.e., the XMM- 
Newton observations). The O Viii intensity predicted 
by our SWCX model is consistent with the SWCX 
O VIII intensity observed by XMM-Newton, while the 
IKoutroumpa et all (|2006l . l2007f) model slightly underpre- 
dicts the O VIII intensity. However, both models under- 
predict the SWCX O vii intensity observed by XMM- 
Newton. 

The implications of these results are discussed in the 
following section. 

7. DISCUSSION 

7.L Solar Wind Charge Exchange in the XMM-Newton 
Spectra - Evidence for a Localized Solar Wind 
Enhancement 

To summarize the preceding sections, we have found 
evidence that the oxygen line intensities in our XMM- 
Newton spectra are enhanced with respect to those in 
our Suzaku spectra. These enhanced intensities are most 
likely due to SWCX emission. Models for heliospheric 
SWCX emission imply that at least part of the differ- 
ence between the XMM-Newton and Suzaku spectra can 
be explain by a change in the global state of the so- 
lar wind between solar maxim um and solar minimum 
(jKoutroumpa et a"ni2006l [20071) . This is because at solar 
maximum there are more of the high-stage oxygen ions 
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which give rise to heliospheric oxygen emission, leading 
to brighter heliospheric oxyge n lines. 

Ho wever, the "ground level" iKoutroumpa et al] (|2006l 
l2007f ) model, which uses typical values for the various rel- 
evant solar wind parameters, underpredicts the SWCX 
oxygen intensities measured with XMM-Newton. This 
implies there is an additional source of SWCX emission, 
in addition to that which is generally expected at solar 
maximum. We suggest that the increased oxygen emis- 
sion during our XMM-Newton observations was due to 
a localized enhancement in the solar wind moving across 
the line of sight. This enhancement could have been 
a region of increas ed density, or a region of enhanced 
oxygen abundance. IKoutroumpa et aLr(|2007D have also 
suggested that there was an enhancement in the solar 
wind during the XMM-Newton observations, which they 
attribute to a coronal mass ejection starting 2.3 days be- 
fore the start of the XMM-Newton observations. 

The "localized solar wind enhancement" scenario we 
have suggested to explain our observations could have 
important implications for determining whether or not a 
SXRB spectrum is likely to be contaminated by SWCX 
emission. Previous observations of SWCX emission in 
the SXRB have found that times of increased SWCX 
emission appear to be associated with enhancements in 
the solar wind as measured by satellites s uch as ACE 
(|Snowden et al.l 120041 : iFujimoto et"an [2007|) . However , 
in the scenario suggested above, the solar wind enhance- 
ment during the XMM-Newton observations was away 
from the Earth, and so not detectable by ACE. As a re- 
sult, the ACE data from around the time of the XMM- 
Newton o bservations show no unusual features (Fig. [8]) , 
which led iHenlev et all (|2007f ) to conclude that SWCX 
was unlikely to be significantly contaminating their spec- 
tra. Indeed, if contemporaneous solar wind data are used 
as an indicator of SWCX contamination, Figure [8] sug- 
gests at first glance that the SWCX emission would be 
higher during the Suzaku observations than during the 
XMM-Newton observations, which is the exact opposite 
of what we observe. Our results therefore imply that 
simply inspecting contemporaneous solar wind data from 
ACE or other satellites might not be sufficient for deter- 
mining whether or not a SXRB spectrum is contaminated 
by SWCX emission. 

7.2. Solar Wind Charge Exchange in the Suzaku 
Spectra ? 

Our analysis suggests that while the XMM-Newton 
spectra are badly contaminated by SWCX emission, par- 
ticularly in the oxygen lines, the oxygen lines in the 
Suzaku spectra are not badly contaminated. Other exist- 
ing XMM-Newton and Suzaku observations yield oxygen 
line int ensities due to SWCX of ^^5-7 L.U. for O v ii and 
O VIII (jSnowden et al.ll2004l; [Fuiimoto et al.ll2007D . The 
upper limits on the SWCX oxygen emission in our Suzaku 
spectra are much less than these values (see Table [5]) . 

However, although the oxygen lines are apparently un- 
contaminated in our Suzaku spectra, there do appear to 
be other lines in the spectra which may be due to SWCX 
(see Fig. [3]). The off-filament Suzaku spectrum shows ex- 
cess emission at ~0.5 keV. For the 6 adjacent channels 
which make up this feature, = 14.6 when compared 
with our preferred model (Model 2 in Table [7]). For 6 
degrees of freedom, this corresponds to a probability 



of 0.024, implying the feature is significant at the 5% 
level. As noted in §3.2) this feature may be too narrow 
to be an emission line. However, if it is an emission line, 
it could be N vil Lya produced by SWCX. The line is 
unlikely to be due to scattering of solar X-rays off at- 
mospheric nitrogen, as the line remains in the spectrum 
even if we exclude times when the elevation of the line 
of sight above the Earth's limb is less than 20° or 25°, 
instead of the default 10° used in fj2] 

The Ne ix feature at ~0.9 keV also exhibits excess 
emission over the best-fitting model - in this case the 
excess emission is seen in both Suzaku spectra. This 
excess emission could be due to an overabundance of neon 
in the halo (see Tabled]). An alternative explanation is 
that the excess Ne ix emission is due to SWCX. 

There is also excess emission in the on-filament Suzaku 
spectrum at ~1.3 keV, which is probably from Mg xi. By 
fitting a Gaussian to this feature, we have measured its 
intensity to be 0.10 ± 0.07 L.U., which implies it is only 
a marginal detection. As mentioned in §3.21 Mg XI emis- 
sion may also be present in the off-filament Suzaku spec- 
trum, where it causes the Gaussian that is modeling the 
Al K instrumental line to shift to a lower energy and in- 
crease in width. The Mg xi emission could originate from 
the halo if the magnesium abundance is enhanced by a 
factor of ~4 relative to oxygen (see Table[2]). However, as 
with neon, an alternative explanation is that the Mg xi 
emission is due to SWCX. SWCX emission from Mg xi 
has been observed i n other XMM-Newton and Suzaku 
spectra of the SXRB. ISnowden et al.l (| 2004h found Mg xi 
emission due to SWCX among a set of XMM-Newton 
spectra of the Hubble Deep Field North, wi t h an inten- 
sity of 0.40±0.08 L.U., while iFujimoto et al.l (|2007t ) mea- 
sured a Mg XI SWCX intensity of Q-lito^Q L.U. from a 
Suzaku spectrum of the North Ecliptic Pole. These other 
results imply that the Mg xi feature in our on-filament 
Suzaku spectrum is not unusually bright for SWCX emis- 
sion. 

The observation of SWCX emission from Mg xi re- 
quires the presence of H-like Mg+^^ ions in the so- 
lar wind. Mg"*"^^ is the dominant ionization stage of 
magnesium in the te mperature range ~7-9 x 10^ K 
(|Mazzotta et al.|[l998D . As the solar wind ions that give 
rise to SWCX originate in the solar corona, and are 
frozen into the plasma as it expands away from the Sun 
(jCravens 2002), the observation of SWCX emission from 
Mg XI may imply very high temperatures in the solar 
corona. However, more detailed modeling is needed to 
determine if the observed Mg xi SWCX intensities im- 
ply an unusually high Mg+^^ ion fraction in the solar 
wind. Such modeling of the Mg xi and other SWCX 
fines is beyond the scope of this paper. Nevertheless, 
the presence of these lines in our Suzaku spectra raises 
an interesting question: if there is SWCX emission from 
Mg XI, and possibly Ne ix and N vii, why is there not 
significant SWCX emission from O vii and O viii? 

7.3. The Normalization of the Extragalactic X-ray 
Background 

As was noted in t j3.11 and as can be seen in Table [71 
the normalization of the extragalactic background is dif- 
ferent in our two Suzaku spectra, with the off-filament 
value (7.7 ±0.3 photons cm^^ s"^ sr"^ keV"^ at 1 keV) 
lower than the expected value of -^10 photons cm^^ s~^ 
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sr^^ keV~^ (e.e.. IChen et al.lll997[ ). One possible rea- 
son for this is incorrect subtraction of the particle back- 
ground: if the particle background spectrum used for 
our off- filament observation is too bright, this would re- 
sult in our background-subtracted SXRB spectrum be- 
ing fainter than expected. However, we do not think this 
is likely to be the cause. At high energies we do not 
expect any cosmic X-ray emission, and so any counts 
in our Suzaku spectra at these energies will be due to 
the particle background. For each observation, the 10- 
12 keV count-rates in the non-background-subtracted 
source spectrum and in the particle background spec- 
trum (generated from night-Earth data) agree within 8%. 
In contrast to this, the particle background flux would 
have to be reduced by ~30% in order to increase the 
off-filament extragalatic normalization to 10.5 photons 
cm^^ s^^ sr^^ keV^^. It therefore does not seem likely 
that an incorrect particle background normalization is 
the reason for this discrepancy in the on- and off-filament 
extragalactic normalizations. 

As well as allowing the extragalactic normalization to 
differ between our two Suzaku spectra, while keeping ev- 
erything else the same, we also tried two other methods 
for dealing with this discrepancy. First, we simply ig- 
nored the discrepancy, and forced the model to be iden- 
tical for both the on- and off-filament datasets. Second, 
we tried multiplying our entire model (Local Bubble -I- 
halo -I- EPL) by a variable factor /og when fitting it 
to our off-filament Suzaku spectrum; this multiplicative 
factor was an additional free parameter in the fit. The 
second method essentially assumes that the discrepancy 
is due to a reduction of ^20% in Suzaku^s sensitivity be- 
tween the on- and off-filament observations. We do not 
think this is plausible, but we still tried out this model. 

We found that these two new methods gave tempera- 
tures in good agreement with each other and with those 
in Table [71 although the emission measures of the Local 
Bubble and of the cool halo component varied by up to 
~15% and ~25% between the models, respectively. Of 
the three methods we tried, our original method (i.e., 
that presented in Table [7]) gave the best fits to the data. 
Also, when wc fitted our multicomponcnt model to our 
on- and off-filament spectra individually, we found that 
the plasma model components are in good agreement 
with each other (Table [71 Models 6 and 7). As was noted 
in this seems to support our method for handling 
the discrepancy in the extragalactic normalization. 

In summary, we do not know why the extragalactic 
normalization differs between our two spectra. It seems 
unlikely to be due to incorrect particle-background sub- 
traction, or to a ^20% reduction in Suzaku^s sensitivity 
in the space of a few days. We have tried several dif- 
ferent methods for handling this discrepancy, and the re- 
sults they give are generally in good agreement with each 
other. We therefore conclude that, despite the discrep- 
ancy with the extragalactic normalization, the plasma 
model parameters we have obtained should be reason- 
ably robust. 

7.4. The Local Bubble 

Because of the SWCX contamination of the oxygen 
lines in the XMM-Newton spectra, and the apparent lack 
of such contamination in the Suzaku spectra, the best- 
fitting model of the LB we have derived from our Suzaku 



spectr a is quite different from that in iHenlev et al.l 
(|2007D . In that project, all the foreground oxygen emis- 
sion in the XMM-Newton spectra was thought to be 
from the LB, leading to a higher LB temperature than 
we have obtained from our Suzaku + ROSAT analysis: 
log(rLB/K) = 6.06, instead of log(TLB/K) = 5.98. The 
lower foreground oxygen intensities in the Suzaku spectra 
cannot be explained by simply lowering the LB emission 
measure, as the ROSAT R12 data place a constraint on 
the LB emission measure. 

Note that our LB emission measure (0.0064 ± 0.0003; 
Table [71) is roughly one-third of the value determined by 
iHenlev et al.l (|2007[ ). However, this is merely because 
we used different plasma emission codes to model the 
ROSAT R12 data. We used t he Raymond- Smith code 
to model the R12 data, whereas fHenlev et al.l (|2007,) used 
APEC. For temperatures around 10^ K, the Raymond- 
Smith code predicts ~3 times as much 1/4- keV flux as 
APEC. This means that for a given observed flux the 
Raymond-Smith-derived emission measure will be ^3 
times smaller than the APEC-derived value, which is ex- 
actly what we find. Our LB emission measure is in rea- 
sonable agreement with the range of LB emission mea- 
sures determined from the ROSA T All-Sky Survey 1/4 - 
keV data (0.0018-0.0058 cm-^ pc: lSnowden et al.H'l998D . 
which was also modeled using the Raymond-Smith code. 

The LB temperature that we measure (log[TLB/K] — 
5.981 

0.04! Table [7]) is lower than the temperatures ob- 
tained from other recent shadowing observations made 
with XMM-Newton and Suzaku: log(rLB/K) « 6.04- 

usin g MBM 20 as a shadow) 
20071 using M BM 12). The 



6.06 (Gal eazzi et al.l 12007 
and «j6.08 ([Smith et al " 



range of values from [Galeazzi are from their 

resul t s for different abundance t ables (I Anders fc Ebihar^ 
i982; Anders & Grevcssc 19891: iGrevesse fc Sauvallll998l : 
Lodders 2003). Their derived LB temperature is fairly 
insensitive to the set of abundances that they used. If 
we repeat our spectral analysis using [Anders fc Grevessd 
( 1989 ) abundances, we obtain log(rLB /K) = 6.07 ± 

0.05, whic h is consist ent with the [Galeazzi et al.l 

([2007( ) and [Smith et all ([2007[ ) temperatures. How- 
eve r, when we us e the other abundances tables used 
by [Galeazzi et al.l ([2007[ ). we obtain LB temperatures 
of log(TLB/ K) = 5.98^^4 ([Anders fc Ebiharal [1982') ■ 
5.97 ± 0. 04 ([Grevesse fc Sauvall [l99l. and 5.96 ± 0.04 
([Lodderd [20031 ). These values are closer to our origi- 
nal LB tempeature, obtained using iWilms et al.l ([200fll) 
abun dances (which were not used bv [Galeazzi et al.l 
[2007[ ). Hence, while the LB tem perature we ob- 
tain using [Anders fc Grevessel ([l989t) abundanc e s is i n 
good agreement with that of [Galeazzi et al.l ([2007D . 
the t emperatures obtained u s ing other abundance ta- 
bles ("Anders fc Ebiharal [1982t [Grevesse fc Sauvall Il998[ : 
Lodders 2003) are systematically lower than those of 
Gal eazzi et al.l (|2007ll . These results imply that the LB 
may not be isothermal. 

If we were to assume that the LB is isothermal, it 
would make it more difficult to explain the lower fore- 
ground O VII intensity measured from our Suzaku spec- 
tra: l.lj: ^ L.U. (see Table [j]) compared with 2.63 ± 
(1.78 L.U. (iGaleazzi et al.1 [200I and 3.53 ± 0.26 L.U. 
([Smith et al.l 120071) . In an isothermal LB, a larger X- 
ray intensity in a given direction can be explained by the 
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LB being of greater extent in that direction. Analysis of 
the ROSAT Ah-Sky Survey suggests that the LB radius 
toward the filament is similar to that toward MBM 20, 
and the LB ra dius toward MBM 12 is roughly two-thirds 
of this value (ISnowden et al.lll998L their Fig. 10). This 
means that the different observed O vii intensities can- 
not be explained by the LB being of different extent in 
different directions. However, if the LB temperature is 
higher toward MBM 20 and MBM 12, this would give 
rise to brighter O vii emission in these directions. 

If we assume the radius of the LB is 100 pc in 
the direction we have analyzed, our emission measure 
gives an electron density of 0.0080 cm~^. Combin- 
ing this with our LB temperature gives a LB ther- 
mal pressure of plb/^ = 1.92 ncTLB = 14,700 cm^^ K 
(13,100-16,100 cm-3 K). Galeazzi et all (120 Q7) obtained 
j PLB/fc = 16,200-19,500 cm-3 K, while ISmith et all 
(|2007f) obtained plb/^ = 22,000 cm'^ K (assuming 
log[TLB/K] — 6.08). These results suggest that the 
LB might not be in pressure equilibrium. Alterna- 
tively, they m a y indi cate that the spec t ra an alyzed by 
iGaleazzi et all (|2007t) and ISmith et al.1 (|2007D are con- 
taminated by SWCX emission. 

7.5. The Halo 

As we found in i )3.3[ a possible explanation for the 
excess Ne ix emission visible at ~0.9 keV in Figures[3]and 
[5] is that neon is enhanced in the halo, relative to oxygen 
and iron. (The alternative is that this emisson is due to 
SWCX, as mentioned in §7.21 ) If neon is really enhanced 
relative to oxygen and iron in the halo (or, equivalently, 
if oxygen and iron are depleted relative to neon), this 
may be evidence that some of the oxygen and iron in the 
halo is in dust, rather than in the hot gas. Neon, being 
chemically inert, does not condense into dust. 

Our best-fitting halo model predicts ~300 L.U. for 
the intrinsic O VI AA1032, 1038 doublet intensity (us- 
ing data from ATOM DB). Note th a t this value is much 
smaller than that in iHenlev et al.l (|2007f ). This is be- 
cause their cooler halo component was cooler than ours, 
and had a larger emission measure (compare Models 2 
and 9 in Table [H). Our model O vi intensity is ~20 
times smaller th an the intrinsic halo O VI int ensity 
determined by Sh elton. Sallmen. fc Jenkini ()2007l ) from 
an off-filament Far Ultraviolet Spectroscopic Explorer 
(FUSE) spectrum. This implies that there is both hot. 
X-ray-emitting gas and warm, ultraviol et-emitting gas in 
the Galactic halo. IShelton et al.l ()2007t ) present a model 
for the distribution of warm and hot gas in the halo, using 
O VI and ROSAT data as constraints. Further modeling 
of the hot gas distribution in the halo, using the Suzaku 
data discussed here, will be presented elsewhere (S. J. 
Lei et al., in preparation). 

8. SUMMARY & CONCLUSIONS 

We have analyzed a pair of Suzaku spectra of the 
SXRB, obtained from pointings toward and to the side 
of a nearby absorbing filament in the southern Galac- 
tic hemisphere. We used a shadowing technique to de- 
compose the spectra into their various thermal emission 
components, due to the LB and the Galactic halo, us- 
ing ROSAT 1/4-keV data to help constrain the model at 
lower energies. We find that our best-fitting model pa- 
rameters (temperatures and emission measures) do not 



agree with those determined bv lHenlev et al.l (|2007f ) from 
a pair of XMM-Newton observations in the same direc- 
tions as our Suzaku observations. We attribute this dis- 
crepancy to SWCX contamination in the XMM-Newton 
spectra, particularly in the oxygen lines, which was not 
taken into account in the original analysis. 

Our results suggest that inspecting contemporaneous 
solar wind data might not be sufficient for determining 
if a SXRB spectrum is contaminated by SWCX emis- 
sion. We have reached this conclusion by comparing our 
observations with models for SWCX emission. The sim- 
ple model of heliospheric and geocoronal SWCX emission 
that we have presented, which assumes an isotropic so- 
lar wind and uses contemporaneous solar wind data from 
ACE and WIND, cannot explain the different amounts of 
SWCX oxygen emission observed by XMM-Newton and 
Suzaku. The more sophisticated heliosp heric emission 
model of iKoutroumpa et al.1 ()2006l . l2007f) . which takes 
into account the change in the global state of the so- 
lar wind between solar maximum, can partly explain 
why more SWCX emission is observed by XMM-Newton 
than by Suzaku. However, this model underpredicts the 
amount of oxygen emission observed by XMM-Newton 
at solar maximum. Instead, the excess SWCX emission 
in the XMM-Newton spectra may have been due to a lo- 
calized enhancement in the solar wind moving across the 
line of sight during the XMM-Newton observations. If 
this enhancement was not near the Earth, it would not 
have shown up in the solar wind data from ACE and 
WIND. 

While it appears that the oxygen lines in our Suzaku 
spectra are not badly contaminated by SWCX emission, 
we cannot be certain of the true level of SWCX emis- 
sion in our spectra. Indeed, it is possible that there is 
SWCX emission from nitrogen, neon, and magnesium in 
our Suzaku spectra. However, if in general the level of 
SWCX contamination is low in our Suzaku spectra, as we 
think is the case, then our best-fitting models of the LB 
and halo should be more accurate than those obtained by 
iHenlev et al.l (|2007f ) . Our best-fitting LB temperature is 
lower than the temperatures obtained fro m other recent 
XMM-Newton and Suza ku observations (IGaleazzi et al.l 
120071 ISmith et alll2007f l. suggesting that the LB is not 
isothermal. Our LB model also yields a lower pressure 
than these other LB models, implying that the LB may 
not be in pressure equilibrium. Our modeling of the halo 
emission, meanwhile, suggests that oxygen and iron may 
be depleted relative to neon, possibly because they are 
partly in a dust phase. 
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TABLE 7 
Spectral Fit Results 



Local Bubble Halo (cool) Halo (hot) EPL norm'' 

logT RM? logT EMy logT KM? On Off A^c'* 

Model Dataset(s)'^ (K) {10"^ cm"'' pc) (K) (IQ-^ cm-*^ pc) (K) (IQ-^ cm~6 pc) Filament Filament (10^® cm'^) x^/dof 



1 


S 


+ R 




7 9+0.4 
'•^-0.3 


6.12 ±0.01 


24.2l^;| 


6.50 ±0.02 


5 7+0.4 
0- '-0.5 


11.1 ± 0.3 


7.6 ±0.3 


0° 


734.24/703 


2 


S 


+ R 




0-^-0.3 


6.11 ±0.01 


34.6129 


o-ou_o 02 


6 5+°-'^ 
"■0-0.5 


11.2 ± 0.3 


7.7 ± 0.3 


0.28 ±0.04 


698.44/702 


3 


s 


+ R 






6.04 ±0.01 


69.8 ±2.7 


6.48 ±0.02 


7 fi+O-'^ 
'■''-o.e 


11.3 ±0.3 


7.6 ±0.3 


0.40 ± 0.04 


1068.63/704 


4 


s 


+ R 


D-UJ_o.o3 


o O+0.3 
».»_0.4 






fi.35±0.01 


12 n+°-8 


ll-"_0.2 


8.6 ±0.3 


< 0.08' 


870.51/704 


5 




S 


5.98° 


1 4+^-2 


6 09+" °2 
D-U9_o.o3 


41 7+3.2 


fi 49+0 °i 


fi 

"■°-0.9 


11.2 ±0.3 


7.6 ±0.3 


0.28° 


672.79/690 


6 


s 


(on) 






fi n9+o °3 

D->JM_o.o4 


4lt« 


fi.52 ±0.02 


7 n+0.6 


10.9 ±0.3 




0.28° 


345.72/374 


7 


s 


(off) 






6.08 ±0.03 


4,4.2+11 


fi.47 ± 0.03 


6 fi+l ° 
"•"-0.8 






0.28° 


317.14/313 


8 




SS 


5.98° 


n "3+5-8 

"■•^-0.3 


fi 0fi+° °3 


47.Stti 


°-^*-0.03 


7 1+1.0 
'■^-0.9 


11.0 ±0.3 


7.6 ±0.3 


0.28° 


655.01/687 


9 


X 


+ 




18 


^ Qo+0.04 


170 


6.43 ±0.02 


11 


10.5° 


10.5° 




435.86/439 


10 


X 


+ R' 


6.30 ±0.01 


12 5+°-® 


p. -0+0.14 
"J- '■^-0.12 


160+3°° 
J^ou_ioo 


fi ^fi+0 06 




11.2° 


7.7° 




442.21/439 



Note. — Model 2 is our preferred model. 
^ Extragalactic power-law normalization at 1 keV in photons cm~^ s~^ sr~^ keV~^, assuming a photon index of 1.46. 
^ Emission measure E.M. = / n'^dl. 
° S = Suzaku; R = ROSAT; X = XMM-Newton. 

Carbon column density of the vphabs model used to model XIS contamination above that included in the CALDB (see i|3.2l for details). 
° Value frozen during fitting. 

' 2a upper limit. 

^ Fit with variable abundances (see i|3.3l and Table [2J. 

From lHenlev et al.l l(200l) . 
' Data reanalyzed to match method used for Suzaku (see 13] for details). 



